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Abstract 



In 1855 H. J. S. Smith [3] proved Fermat's Two Squares using the 
notion of pahndromic continuants. In this paper we extend Smith's ap- 
^^ . proach to proper quadratic form representations in some commutative 

ff^ I Euclidean rings, including rings of integers and rings of polynomials 

C^ ■ over fields of characteristic different from 2. 

While our results are natural extensions of Smith's approach and 
are unlikely to be new, we believe our presentation is new. 
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1 Introduction 

Fermat's Two Squares Theorem, that is. 
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A prime number p is representable by the form x^ + y^ iff — 1 is 
a quadratic residue modulo p. 

is with no doubt a remarkable result. Many proofs of the theorem have been 
provided; see, for instance, [1,3,7,11,15]. It is also true that most proofs have 
much in common, for instance. Smith's proof is very similar to Hermite's [7], 
Serret's [11], and Brillhart's [1]. 

We first recall the definitions of Euclidean rings and continuants. 

Definition 1 ( [8, p. 148]) Euclidean rings are rings R with no zero di- 
visors which are endowed with a Euclidean function N from R to the non- 
negative integers such that, for all mi,m2 G R with mi ^ 0, there exist 
q,r E R such that m2 = qmi + r and N(r) < N(mi) . 

Among well-known examples of Euclidean rings we will make use of the 
integers with N(-u) = \u\ and the polynomials over a field with N(0) = and 
N(P) = 2'^*^srec(P)^ jj^ ^]-^jg po^pQY we ouly consider commutative rings. 

Definition 2 (Continuants in arbitrary rings) Let Q be a sequence of 
elements qi, q2, ■ ■ ■ ,qn of a ring R. We associate with Q an element [Q] of R 
via the following recurrence formula 

[] = 1, ki] = ?i, [Qi, 12] = qiq2 + 1, and 
[gi, g2, • • • , Qn] = [gi, • • • , qn-i\qn + [qi, • • • , gn-2] «/n > 3. 

The value [Q] is called the continuant of the sequence Q. 

Properties of continuants in commutative rings are given in [5, ch. 6.7]. 
Here we restrict ourselves to giving an identity due to Lewis Carroll (alias 
Charles Lutwidge Dodgson). 

Lemma 1 (Lewis Carroll identity) Let C be an n x n matrix in a com- 
mutative ring. Let Cii,...,J3yi,...j^ denote the matrix obtained from C by omit- 
ting the rows ii, . . . ,is and the columns ji, . . . ,js. Then 

det(C) detidj-ij) = det{Ci,i) det{Cj.j) - det(ay) det{Cj,i) 

where det(M) denotes the determinant of a matrix M and the determinant 
of the 0x0 matrix is 1 for convenience. 



The use of Lewis Carroll identity provides more properties, including the 
following. 

P-1 [gi,g2,---,gn][g2, •••,?«-!] = [gi,...,g„-i][g2, ••-,?«] + (-1)" (n > 2). 

Given two elements mi and m2 in an Euclidean ring R, the Euclidean 
algorithm outputs a sequence (gi, 52, • • • , ?n) of quotients and a gcd h of mi 
and m,2- The sequence of quotients given by the Euclidean algorithm is 
called the continuant representation of mi and m2 as we have the equalities 
"^1 = [91; 92, • • • , qn]h and m2 = [g2, • • • , gn]/i unless m2 = 0. 

A representation of an element m by the form Q{x, y) = ax^ + 7x1/ + Py"^ 
is called proper if gcd(a;, y) = 1. In this paper we are mostly concerned with 
proper representations. 

For us the quadratic forms f{x,y) = ax"^ + hxy + cy"^ and g{x,y) = 
Ax"^ + Bxy + Cy"^ are equivalent if there is a 2 x 2 matrix M = (aij) with 
determinant 1 such that g{x,y) = f{aiix + ai2y,a2ix + 0221/)- 

For equivalent forms / and g it follows that an element m is (properly) 
represented by / iff ?72 is (properly) represented by g. 

Let p be a prime number of the form 4k + 1. In his proof of Fermat's 
Two Squares Theorem Smith [3] first shows the existence of a palindromic 
sequence Q = (gi, . . . , g^, g^, . . . , gi) such that p = [Q] through an elegant 
parity argument. This sequence then allows him to derive a solution for 
z^ + 1 = (mod p) and a representation x'^ + y'^ for p. 

Brillhart's optimisation [1] on Smith's construction took full advantage 
of the palindromic structure of the sequence (gi, . . . , g^-i, g^, g^, g^-i, . . . , gi) 
given by the Euclidean algorithm on p and zq. He noted that the Euclidean 
algorithm gives the remainders 

r-i = [qi+2, • • • , qs-i, qs, qs, qs-i, ■■■,qi] (i = 1, . . . , 2s - 1), and 
r2s = 

so, in virtue of Smith's construction, rather than computing the whole se- 
quence we need to obtain 

X = rs-i = [qs,qs-i,---,qi] 

y = Ts = [gs-i,---,gi] 

In this case, we have y < x < y^, Brillhart's stopping criterium. 

Extensions of Brillhart's ideas were presented in [6,12,13], where the au- 
thors provided algorithms, at the expense of the palindromic nature of the 



continuant, for finding proper representations of natural numbers as primi- 
tive, positive-definite, integral and binary quadratic forms. For us the set of 
natural numbers N includes the zero. 

This paper can be considered as a follow-up of our earlier paper [4]. In [4] 
we studied the use of continuants in some integer representations (e.g. sums 
of four squares) and sums of two squares in rings of polynomials over fields 
of characteristic different from 2. Here we deal with proper representations 
x^ -|- gxy + hy'^ in Euclidean rings via modified continuants. In Section 2 
we modify the notion of continuant and describe some of their properties. 
Section 3 is devoted to studying proper representations x'^+gxy+hy^ in some 
commutative Euclidean rings, mainly in the ring of polynomials over a field of 
characteristic other than 2. In Section 4 we consider proper representations 
x^ -|- gxy + hy'^ in the ring of integers. Some final remarks are presented in 
Section 5. 

While the results presented here are unlikely to be new, we believe our 
presentation and extension of Smith's ideas are new. 



2 Modified Continuants 

With the aim of considering the problem of properly representing an element 
m as x"^ + gxy + hy'^, we extend the notion of continuants. 

Definition 3 (Modified Continuants in Arbitrary Rings) In a ring R 
associate with the element [Q; h, s] the 3-tuple formed from a sequence Q of 
elements qi,q2, . ■ . ,qn of R, an element h of R and an integer s via the 
following recurrence formula 

{[qi,...qn] ifs>n 

[gi, . . . g„-i]g„ + [gi, . . . qn-2]h if s = n-l 

[gi, . . . g„_i; /i, s]g„ + [gi, . . . g„_2; h,s\ if s < n - I 

The definition of modified continuant carries several consequences, which 
we will refer to as Modified Continuant Properties. 

P-2 [gi, . . . , g„; /i, s] = [gi, . . . , g,_i]/i[g,+2, . . . , g„] + [gi, . . . , qs][qs+i, • • • , Qn]- 
Divide the products of subsequences oiQ = (gi, g2, . . . , g™) obtained by 
removing disjoint pairs of consecutive elements of Q into two groups, 
depending on whether or not the products remove the pair qsqs+i (1 < 
s <n). 



u 



'^qi,uq2,...,u^ ^)"g„;/i,s] 



P-3 If in a ring R we find a unit u commuting with all gj's, then 

[gi, . . . , g„; h, s] for even n 

u~^[qi, . . . ,qn',h, s] for odd n 

P-4 Induction can provide us with the following. 

[gi , . . . g„] if s > n or s = 

qi[q2, ■ ■ ■ g„-i, qn] + h[q3, . . . g„] if s = 1 

[gi, . . . qn-i]qn + [gi, . . . qn-2]h if s = n - 1 

qi[q2, ■ ■ ■ qn] h, s - 1] ii2<s<n-l 

^ + [g3,.-.g„-2;/i,s-2] 



[qi,...qn;h,s] = < 



2.1 Modified continuants in commutative rings 

P-5 The modified continuant [gi, . . . ,qn',h, s] is the determinant of the tridi- 
agonal n x n matrix A = (ajj) with Oj^j = g^ for 1 < i < n, aj^j+i = 1 
for 1 < i < n, tts+i^s = —h and flj+i^j = — 1 for 1 < i < n and % ^ s. 
See the determinant of the matrix below for a small example. 



[gi,g2,g3,g4,g5;/i,3] = det( 



91 1 

-1 92 1 

-1 g3 1 

— /i g4 1 
-1 gs 



P-6 [gi, g2, . . . , g„; /i, n- s] = [g„, . . . , g2, gi; /i, s]. This follows from applying 
Property P-2 on both sides of the equality. 



3 From Q(x^ y) to Q{z^ 1) and back 

In this section, considering the form Q{x, y) = x^ + gxy + hy'^, we deal with 
the problem of going from a representation Q{x,y) of an element m to a 
multiple Q{z, 1) of m and back. We begin with a very general remark which 
is valid for every commutative ring. 

Proposition 1 // Rx + Ry = R then there exist z & R such that Q{x,y) 
divides Q{z, 1), where Rm denotes the ideal generated by m. 

If R is Euclidean, we can explicitly find z and the quotient Q{z, 1)/Q{x, y) 
with modified continuants. 



Proof. We have u and v such that xu + yv = 1. Then, computation 
with norms in the ring obtained from R by adjoining formally a root of the 
polynomial T^ — gT + h provides the identity 

Q{x, y)Q{y — ug, u) = Q{xv — xug — yuh, xu + yv), 

which proves the first assertion. This identity can be interpreted also as a 
kind of Lewis-Carroll identity. 

The determinant of the tridiagonal matrix 



Qn 



12 



M 



-1 qi ••• 
-h qi+g 
-1 



?2 



-1 Qn 

. , qn] and y = [q2, . . . , g„] if n > 1. 

g„_i], [gs, • • • qn-i]) = Q{z, 1), where z = (-l)"+ic 
and c is the determinant of the matrix formed by the 2n — 1 first rows and 
columns of M. D 



is Q{x,y) with x = [qi, 
Also, Qix,y)Q{[qi, 



Remark 1 Let R be a commutative ring. 

If 2 is invertible, the form x^ -\- gxy + hy"^ can be rewritten as {x + gy/2)'^ + 
{h — g^/4)y'^. We may then assume g = without loss of generality. 

If moreover —h is an invertible square, say h + k"^ = 0, then x = (^^) + 

Below we provide a proposition which can be considered as an extension 
of [4, Proposition 2]. 

Proposition 2 Let R = ¥[X] be the ring of polynomials over a field F with 
characteristic different of 2, and —h a (non-null) non-square of¥. 

If m divides z^ + ht^ with z, t coprime, then m is an associate of some 
x^ + hy^ with x, y coprime. 
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Proof. We introduce the extension G of F by a square root u of —h. The 
ring G[X] is principal and z^ + ht^ factorises as {z~ujt){z + ut), with the two 
factors being coprime. Introduce x + uy = gcd(m, z + ut); then x — uy is a 
gcd of m and z — ut, using the natural automorphism of G. The polynomials 
X — uy and x + uy are coprime and both divide m. Thus, m is divisible by 
{x — uy) {x + uy) = x^ + hy'^ . On the other hand, m divides {z — ut){z + ut) . 
Consequently, m is an associate of x"^ + hy'^. Since x — uy and x + uy are 
coprime, x and y are coprime. D 

Proposition 3 Let R = ¥[X] be the ring of polynomials over a field ¥ with 
characteristic different of 2, and let h be a polynomial of degree 1. 

If m divides z^ + ht^ with z, t coprime, then m is an associate of some 
x^ + hy^ with x, y coprime. 

Proof. Consider the extension of the ring R = ¥[X] by a root of T^ + h; 
this extension of R is isomorphic to ¥[T]. 

If m is a multiple of h, then m is an associate of (hy)'^ + hx'^, with hy and 
X coprime. This representation of m follows from a representation x^ + hy"^ 
of m/h. 

We may now assume that h and z are coprime. Then, z'^ + ht"^ factors 
as {z — Tt){z + Tt), with the two factors being coprime. Reasoning as in 
Proposition 2, we let x + Ty be the gcd of m and z + Tt and we obtain that 
X — Ty is the gcd of m and z — Tt and that m is an associate of x^ + hy"^ 
with X, y coprime. D 

Propositions 2 and 3 give some cases where a reciprocal of the first assertion 
in Proposition 1 holds. Next we show how modified continuants and these 
propositions can be used to find the representations. 

The next remark generalises [4, Remark 3]. 

Remark 2 (Algorithmic considerations) For the cases covered in Pro- 
positions 2 and 3 we can resort to Brillhart's [1] optimisation to accelerate 
the computation. We stop the division of m by z when a remainder Ts-i with 
degree at most deg(m)/2 is encountered. This will be the {s — l)-th remainder, 
and {uqs,u~^qs-i, . ■ . ,u^^^^" ^2) o'^e ihe quotients so far obtained. Then 



X 



Ts-i for odd s 

M^Vs_i for even s 



[uqs, u ^qs-i, ...,«( ^)° '^2] for odd s 

u'^uqs, u~^ Qs-i, ... lu''^^^' ^2] for even s 

This remark follows from dividing m/u = [g^, . . . ,qi,qi, . . . ,qs;h,s] by z = 
[qs-i, . . . , qi, qi, . . . , qs] h, s — 1] using modified continuant properties. 

The argument brought forward in [4, Proposition 3] can be apphed to the 
form x"^ + hy"^ in polynomials over a field F of characteristic different from 2, 
where /i G F is a non- square. 

Corollary 1 (of Proposition 2: h non-square unit in F[X]) Let m be 

a non-unit of¥[X] and a divisor of z^ + h for some z G F[X] with deg(2;) < 
deg(m). Then, m = {x^ + hy'^)u for some unit u and the Euclidean algorithm 
on m and z gives the unit u and the sequence 

{uqs, W^qs-i, . . . , M^^^^'^'gi, M^^^^'/i^^gi, . . . , u'^h^'^'^° q^) 

such that X = [qi, . . .qs] and y = [q2, . . . ,qs]. 



In the next example we illustrate Remark 2 and the method of Corollary 1, 
in this order. Let h = 3 and m = 1 + 2X + 3X^ + 2X^ + X'^, then m divides 
((5 + 12X + QX^ + 4X3)/3)2 + 3. The Euchdean division gives 

1 + 2X + 3X2 + 2^3 + x^ =((5 + 12X + QX^ + 4X3)/3)(3X/4 + 3/8) 

+ 3/8 - 3X/4 - 3x74 

Here the first remainder has degree at most deg(m)/2, thus we stop the 
division process and obtain s = 2, x = (3/8 — 3X/4 — 3X'^/4:)/u and y = 
[3X/4 + 3/8] /m. It is now plain to get u = 9/16. 

If instead we use the method of Corollary 1, then we obtain the unit 
u = 9/16 and the sequence (9/16 ■ 2/3 ■ (1 + 2X), 16/9 ■ (-1/2 - X), 9/16 • 
1/3 ■ (-1/2 - X), 16/9 ■ 3 ■ 2/3 ■ (1 + 2X)). From this sequence we know that 

x = [2/3-(l + 2X),-l/2-X] 
y = [2/3-{l + 2X)] 

Corollary 2 (of Proposition 3: h of degree 1 in F[X]) Let m be a po- 
lynomial over F[X] and a divisor of z^ + h for some z G F[X] with deg(z) < 
deg(m) and z, h coprime. Then, m = {x^ + hy^)u for some unit u and the 
values of x and y can be obtained by Remark 2. 



Consider the following example, leth = X,m = l + X + X^ + X'^ and 
z = {X^ + 2X2 _^ iy2. Then, the division gives 

1 + X + X^ + X^ =(X3 + 2X2 ^ ^y2 ■ (-2 + 2X) + 2 + 2X2 

At this step we should stop the division process as the first remainder has 
at most half the degree of m. Now we know that s = 2, x = (2 + 2X'^)/u 
and y = u~^[—2 + 2X] for a unit u. It plainly follows that -u = 4. 

We now may wonder how far can we push this method for polynomials 
over a field of characteristic different from 2? That is, will the method work 
for h with deg(/i) > 1 over any such field? 

We first note that the property 

"m|z2 -\- h ^ 3x, y,u {m = u{x^ + y'^h) A u unit)" 

does not hold in general for h reducible. Consider h = X^ + X^ + X in 
polynomials over a field of characteristic 7^ 3. Then, X^ + X + 1 divides 
0^ + l^/i and is certainly not of the form x^ + y'^h. Indeed, here we have either 
y'^h null or of odd degree > 3. In the former case, it follows that x'^+y'^h = x^ 
is a square, but X^ + X + 1 is not a square in a field of characteristic 7^ 3, 
while in the latter case x^ + y'^h has degree > 3 > degree (X^ + X + 1). 

What about irreducible h with deg{h) > 2? Well, already for degree 2 the 
property does not hold in general. Indeed, consider in Q[X] the polynomials 

h = X^-2, z = X^ andm = X-l. Observe that X - 1 divides X^ + X^ - 2 
and that X — 1 is not of the form ^(x^ + y^{X^ "2)). To see this note that 
the degree of x^ + y^{X^ — 2) is either 2 deg(x) or 2 + 2 deg{y). 

For specific fields we find situations where the property holds. Take, for 
instance, the field M of reals, h = X^ + 1 and every real polynomial m taking 
only positive values over R. It is known that the polynomials m over M that 
take at every point of M a positive value has the form fKa^X^ + 2bkX + Cfc), 
where ak,bk,Ck G M and b^ — akCk < 0. Thus, it suffices to consider the case 
of m = aX2 + 26X + c with a > 0, c> and b'^ ~ ac <0. If b = 0, then 

'(v/a^cX)2 + y3^(X2 + l) if a>c 
^/c-d^ + ^/a^{X^ + 1) if a < c. 



If instead & 7^ 0, then, setting d = J {a + c — ^J{a — c)^ + Ab'^)/2, we obtain 

m = {y/a — d'^X + ey/c — d'^y + d'^{X'^ + 1), where e = ±1 has the sign of b. 

For h = — X^ — 1 and every real polynomial m over M, we have another 
situation where the property holds. Observe that the form Q{x,y) = x'^ + 



(— X^ — l)y''^ is equivalent to the form Q{x, y) = x^ + 2Xxy — y"^ (by Remark 
1). Any polynomial of degree 1 is an associate of some a^ — 6^ + 2abX with 
units a and b. We now take care of polynomials m = k{X'^ + 2vX + w) with 
no real zeros and k,v,w E M. Here note that v^ < w. Set p{X) = (X + a)^ + 
26(X + a)X — 6^. We solve the equation p{X) = m in (a, b, k). We first find 
that A; = l + 2b, a = kv/{l + b) (if 6^ -l)andw; = {a^-b'^)/k = -X^-2vX. 
If 6 = — 1 then v = 0, k = —1 and a = ±-\/l — w with < w < 1. If instead 
6 7^ —1, then, substituting a = kv/{l + b) into —(1 + 6)^p(X), we obtain 

6^ + 2(w + 1)6^ + (5w - 4t;^ + 1)6^ + 4(w - t;^)^ + w-v'^ = Q. 

This equation in 6 is 1/16 when b = —1/2 and — f^ when 6 = — 1. Hence 
there is a solution b in the open interval (—1, —1/2) for w > 1. Consequently, 
each real polynomial is an associate of some polynomial x'^ + 2xyX — y^ = 
(x + Xyy + i-1 - X^)y\ 

Using the automorphisms of M.[x], both previous approaches can easily be 
applied to any real polynomial of degree 2 with no real roots. 



3.1 From Q{z,l) to Q{x,y) with assumptions 

Consider the form Q{x, y) = x'^ + gxy + hy"^, now we deal with the following 
question: 

Given z,m in a commutative ring R such that m = Q[x,y)u 
(with u a unit) is a divisor of Q{z,l), under which conditions 
can we recover a representation Q{x,y)u of m from m and z via 
modified continuants? 

Throughout this section we suppose that the element m admits a proper 
representation Q{x,y)u with u a unit. 

By Proposition 1 there exists a multiple of m of the form Q{z,l) with 
z = {—l)^^^[qs_i, . . . , gi +g, qi, . . . , Qs] h,s — 1]. We may then assume m/u = 
[qs, ...,qi+g,qi,...,qs;h, s] with x = [gi, . . . , g^] and y = [q2, ...,qs]- Ob- 
serve that in the case of s = we deal with the form Q{x,y) = x^ + y"^, 
where m/u = [qs, . . . , gi, gi, . . . , g^] and z = [q^-i, . . . , gi, gi, . . . , g^]; this case 
was studied in [4]. 

Conditions for the aforementioned question will heavily depend on the 
ring under consideration. In a ring R such conditions will provide determin- 
istic algorithms for solving the equations m = Q{x,y)u. 
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One way of attacking this problem is by ensuring the uniqueness of the 
division process of m by z. In this context the conditions we are after are 
those guaranteeing that the division process oi m hj z works up to we find 
an associate of x as a remainder. By properties of modified continuants, we 
know that an associate of x crops up as one of the remainders r^-i in the 
division of m by z. Hence, a suitable algorithm goes as follows: divide m by 
z and stop the division at "the right remainder r^.i". If we succeed to do 
that, then 



X 



rs-i 



for odd s 



u ^Tg-i for even s 



[uqs,u ^qs^i,..., 



-u 



-uq. 



u 



(~iy 



m 



for odd s 



s,-u ^qs-i,. . .,-u^ ^^° ^2] for even s 



Here {—iy~^^uqs, {—ly^^u ^qs-i, • • • , {—ly^^u^ ^^^ q2 are the quotients 
so far obtained. This is a generalisation of Remark 2. 

Consequently, we need a criterium to stop the division at the right re- 
mainder Tg^i. 

As an illustration, consider the example below for s = 4. Here m can 
take the form [^4, ^3, q2, qi + g, gi, g2? Q's, Q'4; ^, 4] and there is a 2; of the form 
"bs) Q'25 Q'l + fi'; Q'l; Q'25 93, Q'4; ^5 3]. Then, dividing m by z, in virtue of the 
assumed uniqueness of the division process, we find 

[94, 93, 92, 91 + g, 91, 92, 93, 94; h,A]=- [gg, gs, 9i + ^, 9i, 92, 93, 94; h, 3] ■ (-94^) 

+ [92, 9i + ^, 91,92,93, 94; /i,2]M 
-[93, 92, 9i + 9, 9i, 92, 93, 94; h, 3] =[^2, 9i + ^, 9i, 92, 93, 94; h, 2\u ■ {-qs/u) 

- [9i+^,9i,92,93,94;/i, 1] 
[92, 9i + 9, 9i, 92, 93, 94; h, 2]u = - [q^ + g, q^, q^, q^, q^, h, 1] ■ (-92^) 

+ [91,92,93,94]-" 

Here we should stop as [gi, g2, 93, 94]^ is an associate of x. 

Thus, we first need conditions guaranteeing that the above division pro- 
cess is genuine, which means that 

N(-[93, 92, 9i + 9, 9i, 92, 93, 94; h, 3]) >N([g2, 9i + 9, 9i, 92, 93, 94; h, 2]u) > 
N(-[9i + 9, 91, 92, 93, 94; h, 1]) >N([gi, g2, 93, 94]m)- 
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and secondly, we need a criterium telling us where to stop the division pro- 
cess. 

We consider the case of the polynomials over F[X] with F being of char- 
acteristic other than 2. By Remark 1, we assume Q{x,y) = x'^ + hy^. The 
conditions here are the following. 

(i) Uniqueness of the division process. The Euclidean function N(0) = 
and N(P) = 2'i^s''^<=(^) for P e ¥[X] does guarantee this. 

(ii) The chosen z satisfies N(z) < N{m). This implies that deg(gj) > 1 for 
i = 2, . . . ,s. There is no condition on qi. 

(iii) h is square-free. 

Using properties of modified continuants it is not difficult to see that the 
above division process is genuine. 

We now turn our attention to the stopping criterium. The stopping cri- 
terium is < degree(?72)/2, that is, we stop when we find a remainder Vg-i 
of degree < degree(?72)/2. We now prove this assertion, but first we remark 
that, if /i is a non-square unit in F[X], then we may assume deg(x) > deg{y). 

In the division process we find that the remainder rs-2 is an associate of 
qix + hy and that deg(rs_2) > deg(m)/2. Indeed, if deg(a;) > deg{y) then 
deg(gi) > 1 and we are done. If instead deg{y) > deg(a;) then qi = and 
we only run into "trouble" if deg(/i) = 0, which is actually fine by the above 
remark. It only remains the case of deg(x) = deg{y), but again, we only run 
into "trouble" if deg{h) = 0. Thus, we have proved the assertion. 

A more complex case is Z. For the form Q{x, y) = fx"^ + hy"^ over i? = Z, 
the conditions 

(i) Uniqueness of the division process. To ensure this we restrict the divi- 
sion to N, 

(ii) f,h>l, gcd(/, h) = 1, and m > f + h + 1 with gcd(m, fh) = 1, and 

(iii) 1 < z < m/2. 

together the stopping criterium Vg-i < y^m/f provide the desired elements; 

see [6]. 
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4 From Q{z,l) to Q{x,y): integral quadratic 
forms 

In Subsection 3.1, we presented a deterministic algorithm that, under certain 
conditions, can find representations Q{x, y) for an associate of an element m. 
Such an algorithm assumed the representations Q{x,y) existed. In this sec- 
tion, we do not make such assumptions. Given integers m, z such m\Q{z, 1), 
we provide an algorithm that proves the existence of representations Q{x, y)u 
of m for a unit u and certain forms Q. 

Since 2 is not invertible in Z, we have to consider the rings of algebraic 
integers of Q[\/—h], that is, the rings Z[^/— /i] for forms x^ + hy"^ with \h\ 
square-free and h ^ —1 (mod 4), and the ring Z[(l + ^/l — 4/i)/2] for forms 
x'^ + xy + hy'^ with |1 — Ah\ square-free. 

What are those rings of integers for which the following property holds? 

''m\Q{z, 1) =^ 3x, y,u {m = uQ{x, y) Au unit)." 

The answer is given by the rings whose corresponding forms have class 
number H{A) equal to 1 [2, pp. 6-7,81-84]. In the case of A < 0, all the 
principal rings satisfy the property. These A's are called Heegner's numbers 
(1, 2, 3, 7, 11, 19, 43, 67, 163). In this context, modified continuants provide a 
constructive proof of the property. For the case of A > 0, while we do not 
even know whether the list of such determinants is infinite, it is conjectured 
this is likely the case [2, pp. 81-82]. 

The class number H (A) gives the number of equivalence classes of integral 
binary quadratic forms with discriminant A. 

Below we present a division algorithm (Algorithm 1) which, for any Heeg- 
ner's number d, gives a proper representation Q{x,y) of m assuming m di- 
vides Q{z, 1). Unsatisfactorily, our algorithm does not terminate correctly 
for all A > with H{A) = 1. 

Since m = 1 trivially admits proper representation (1,0) of Q{x,y) = 
x"^ + gxy + hy'^, Algorithm 1 assumes \m\ > 1. 

Remark 3 (Algorithm 1: Prioritising non-null quotients) In the Eu- 
clidean division of Zg-i by rUs with \zs-i\ < \Tns\, a valid quotient kg could 
be ±1 or 0. By "prioritising non-null quotients kg" we mean that, in this 
situation, we always choose the non-null kg- 

Proposition 4 (Algorithm 1 Correctness) Let h,d,u,mQ,Zo andmi,Zi, 
Qi (i = 1, . . . ,s) be as in Algorithm 1. Then, Algorithm 1 produces a proper 
representation Q{x,y) ofmo/u. 
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Algorithm 1: Deterministic algorithm for constructing a proper rep- 
resentation Q{x, y) = x^ + gxy + hy'^ of an element m 

input : A Heegner's number d, an integer rriQ with 1 < |mo| and a 
solution Zq of Q{z, 1) = (mod itiq) with 1 < \zq\ < |mo|. 

output: A proper representation Q{x,y) of itlq/u {u = ±1). 

while \ms\ 7^ 1 do 

s ^ s + 1; 

rus ^Q{zs^i,l)/'ms-u 

find kg, Zg E R such that Zg^i = kgrng + Zg with \zs\ < l^sh 

/* We prioritise non-null quotients ks- */ 

end 

/* Here we have the unit nig and sequence {ki,...,ks). */ 

/* To keep consistency with the previous sections of the 

paper we reverse the subscripts of the quotients */ 
ill, 12, ■ ■ ■ , Qs) ^ {ks, ks^i, . . . , ki); 

return (x, y) 

Proof, li d = —1 (mod 4), we consider the form Q{x,y) = x'^ + xy + 
hy'^ with h = {1 + d)/4, otherwise the form Q{x,y) = x"^ + hy'^ with h = 
d. As the proof method is the same in both cases we restrict ourselves 
to the former case, that is, to the case oi d = 3,7,11,19,43,67,163 and 
h = 1,2,3,5,11,17,41. 

Claim 1. Algorithm 1 terminates with the last rrij being ±1. 

We first prove that Algorithm 1 terminates. As a general approach we 
show |mj+i| < \mi\ (i = 0, . . . , s — 1). Once this decreasing character fails we 
show that the algorithm anyway stops. 

Recall we have \mi\ > \zi\ + 1 for ? = 1, . . . , s — 1. 

Case d = 3,7: \mi\\mi\ > zf + 2\zi\ + 1 > \zf + zi + h\ = |mj||mj_|_i|, 
and thus \mi\ > \mi+i\ for |zj| > 1. In the case that d = 7 and a certain 
Zi, say Zs-i, is equal to 1, we have that 1 + 1 + 2 = nis^ims. From this we 
get that either |ms_i| = 2 and \ms\ = 2 or |ms_i| = 4 and \ms\ = 1. The 
configuration |ms_i| = 4 and \ms\ = 1 will cause the algorithm to stop with 
rUs being a unit. In the case of Im^.il = 2 and \ms\ = 2, in the next step 
we have Zs-i = rrisks + Zs and the algorithm would obtain Zg = —1, which 
implies |ms+i| = 1. 
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Consequently, in these two cases Algorithm 1 terminates with the last rrij 
being a unit. 

Case d = 11, 19,43,67, 163: Analogously, we have that |mj||mj| > zf + 
2\zi\ + 1 > \zf + Zi + h\ = |mj||mj+i|, unless — [^^^J < Zi < h — 1. The 
case of — [^^^J < Zg-i < /i — 2 implies Im^-il > \ms\ with \ms\ = 1 (as 
Q{zs-i, 1) = -2^-1 + ^s-i + /i is a prime). The fact that Q{zs-i, 1) is prime 
follows from a well-known theorem by Rabinowitsch [10]: For a discriminant 
A = 1 — 4k, < —7, H{A) = 1 iff x^ + X + K attains only prime values for 
<x< /t-2. 

// instead Zg^i = h — 1, then \ms-i\ = \ms\ = h, which implies that 
Zs = —1 and Im^^^il = 1, causing the algorithm to stop. 

Now we prove that the last rrij must be ±1. At some point we obtain 
Zs-i = rUsks with Q{zs-i, 1) = rris-ims and nig a divisor of h. Suppose rus 
is not a unit, then it follows that rUg = h (as h is prime), Zg = 0, and that 
rrLg^i is necessarily a unit. 

Claim 2. Algorithm 1 produces a proper representation Q{x,y). 

It only remains to prove that x and y have the claimed expression. First 
we reverse the subscripts of the quotients, that is, the quotient kg becomes qi, 
the quotient kg_i becomes q2, and so on. Thus, after the while loop we have 
that Zg^i = nigqi, where rrig is a unit. We know that Q{mgqi, 1) = rrig^img. 
Consequently, rrig^img = [rrigqi + l,mgqi;h,l]. Then, by Property P-6 and 
Property P-2 it follows 

Zg^2 = [m^^q2,mgqi + l,mgqi;h,2]. 

Then, from the equation mg_2'n^s-i = Q{zg-2-, 1) we obtain 

mg_2mj^ = [m^^q2,mgqi + l,mgqi,m~^q2;h,2] = Q{[mgqi,mJ^q2],[mJ^q2]). 

Continuing this process, we have 

Zq = [m^"^)' \g, K, rngQi + 1, m^gi, K~^; h, s] 
mQm[~^^' =[m[~^^' qg,K,mgqi + l,mgqi,K~^,m[~^'>'' qg;h,s] 

where K = rrig qg^i, . . . ,m~q2 and K~ = m~q2, . . . ,mg qg-i. 

Consequently, from Property P-3 it follows that monig = Q(x,y), 



where 



il > qg^i,ml > qg 

r -1 (-1)""^ (-1)^"^ 1 
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As a result, each root Zq of Q{z, 1) = (mod mo) with 1 < \zq\ < |mo| 
gives rise to a proper representation of moms as Q(x,y). The coprimal- 

ity of X and y follows from Property P-1. D 

Let us see now an example. For the ring Z[(l + v^— 19)/2] the form is x^ + 
xy + 5y^. Noticing 251 ■ 11 = 52^ + 52 + 5, the division gives 



251- 11 

11-7 

7-1 



52^ + 52 + 5 
8^ + 8 + 5 
1' + 1 + 5 



-^ 
-^ 
-^ 



52 



4-11 + ^ 
1-7+1 
1- 1. 



Thus, we have m^ = 1 and {q3,q2,<li) = (4, 1, 1). From this we recover the 
continuant representation of mo = 251 



4 1 
-1 1 1 

-1 1 1 

-5 1 + 1 1 

-1 1 1 
-1 4 



251 = det( 



concluding that 251 = x"^ + xy + 5y'^ with x = [1, 1,4] = 9 and y = [1,4] = 5. 



5 Final remarks 



In Algorithm 1 we require m^ to be ±1. However, this may be an unneces- 
sarily strong restriction. If in Algorithm 1 we replace the condition of the 
while loop by Zg 7^ 0, then this modified Algorithm 1 may also end with 
the last rrij, say nig, being different from ±1. Further, if such nig admits a 
representation as Q{x,y), then the formula 



(x^ + gxy + hy'^){z'^ + gzw + hw"^) ={xz — hywY + 9 



xz 



hyw) 



X 



X {xw + yz + gyw) + 
+ h{xw + yz + gywY 

will provide a desired (possibly rational) representation of m = ttiq for a 
larger number of forms Q{x,y). First recall that in Algorithm 1 



mom 



i-iy 



-x^ + gxy + hy^ . 
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where 

r -1 (-1)'"^ (-1)'"-^ 1 

y=[m^ q2,...,m\ > qs^i^m\ (ls\ 

Then, to recover the representation of mo (associated with zq) we just 
need to expressed mi as Q{x,y). 

This simple modification of Algorithm 1 will provide proper representa- 
tions Q{x, y) of ±m for some forms Q with discriminant A > and -ff (A) = 1 
(see [2, pp. 6-7,81-84]). Below we show some of such forms in the format 
(Q, {list of possible ms's}). Note that iris is a divisor of h. 



\x 



2 



2y\ {±1, ±2}) (x2 + xy- Ay\ {±1, ±2, ±4}) 

[x^ - ?,y\ {±1, ±3}) (a;2 + xy - ly\ {±1, ±7}) 

[x^ - Qy^, {±1, ±2, ±3, ±6}) {x^ + xy - V, {±1, ±3, ±9}) 



[X 



2 



7y^ {±1, ±7}) (x2 + xy- lOy^ {±1, ±2, ±5, ±10}) 



(a;2 + xy- y^, {±1}) (a:^ + xy - ISy'^, {±1, ±13}) 

(a;2 + xy- 3y^, {±1}) (a;^ + xy - 15y^, {±1, ±3, ±5, ±15}) 

If rus is not representable as Q{x,y), then the modified Algorithm 

1 would provide a proper representation of a (possibly rational) multiple of 
m; this is the case of Q{x,y) = x"^ — 6y^, mo = —37410 and zq = 1326, 
where Algorithm 1 would give m^ = 2, s = 2 x = 141 and y = 14. Thus, 
37410 ■ 1/2 = 1412 _ g . ^42^ 

From the previous remarks we have that the modified Algorithm 1 covers 
the cases studied in [9]. Observe that the form x^ — by"^, studied in [9] and 
associated with the non-principal ring Z[v^], has been superseded by the 
form x"^ + xy — y^ associated with the integral closure of Z[v5], that is, 
Z[(l±V5)/2]. 

Furthermore, a minor modification of Algorithm 1 can also provide alter- 
native proofs to Propositions 2 and 3. 

The approach presented in the paper is likely to work for other represen- 
tations if new modified continuants are defined. 

Mathematica® [14] implementations of most of the algorithms presented 
in the paper and other related algorithms are available at 

http : //guillermo . com . au/wiki/List_of .Publications 

under the name of this paper. 
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